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ABSTRACT 

We examine whether active galaxies obey the same relation between black hole mass and stellar velocity dis¬ 
persion as inactive systems, using the largest published sample of velocity dispersions for active nuclei to date. 
The combination of 56 original measurements with objects from the literature not only increases the sample from 
the 15 considered previously to 88 objects, but allows us to cover an unprecedented range in both stellar velocity 
dispersion (30-268 km s -1 ) and black hole mass (10 5 —10 8 6 M 0 ). In the Mbh - cr* relation of active galaxies we 
find a lower zeropoint than the best-fit relation of Tremaine et al. (2002) for inactive galaxies, and an intrinsic 
scatter of 0.4 dex. There is also evidence for a flatter slope at low black hole masses. We discuss potential con¬ 
tributors to the observed offsets, including variations in the geometry of the broad-line region, evolution in the 
Mbh — cr* relation, and differential growth between black holes and galaxy bulges. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: Seyfert 


1. INTRODUCTION 

Black holes (BHs) are a basic component of galaxies, and the 
existence of a tight correlation between the stellar velocity dis¬ 
persion of bulges (cr*) and the BH mass (the Mbh-ct* relation; 
Gebhardt et al. 2000a; Ferrarese & Merritt 2000) suggests that 
the growth of the BH plays a fundamental role in the growth 
of the bulge, although exactly how remains unclear (e.g.. Silk 
& Rees 1998; Kauffmann & Haehnelt 2000; Di Matteo et al. 
2005). The Mbh-ct* relation for inactive galaxies, in so far as 
it represents the final state of the BH-bulge system, represents 
a boundary condition for various evolutionary scenarios, and 
some clues are embedded in the scatter and possible skewness 
in the local Mbh-ct* relation (Robertson et al. 2005). Unfortu¬ 
nately, the number of available points in the relation for inactive 
galaxies is limited, and statistics are poor. Currently, our only 
recourse is to rely on BH masses from active galactic nuclei 
(AGNs). To the extent that the BHs in AGNs continue to gain 
mass, the relation of BH masses in AGNs to their host bulges 
may carry additional information about the establishment of the 
relation for inactive sources. Since the majority of BH mass 
was assembled at high redshift (e.g., Yu & Tremaine 2002), we 
might expect to find the strongest evidence for evolution in the 
Mbh-ct* relation at large cosmological distances (Shields et al. 
2003; Treu et al. 2004; Walter et al. 2004; Peng et al. 2005). 
While undoubtedly this is a vital direction to pursue, there are a 
number of compelling reasons to study the local AGN Mbh - cr* 
relation as well. 

For one thing, while various methods are available to char¬ 
acterize the bulge potential, virial mass estimation (e.g.. Ho 
1999; Wandel et al. 1999; Kaspi et al. 2000), where the dense 
broad-line region (BLR) gas is assumed to be on Keplerian or¬ 
bits around the central BH, is currently the most widely uti¬ 
lized tracer of BH mass. In the absence of a detailed model 
of the BLR, the zeropoint for the virial BH mass scale is set 
through a direct comparison with stellar velocity dispersions 


for a small sample of local AGNs (Gebhardt et al. 2000b; Fer¬ 
rarese et al. 2001; Nelson et al. 2004; Onken et al. 2004). The 
virial BH masses considered in these studies are remarkably 
consistent with the Mbh - cr* relation of inactive galaxies, sug¬ 
gesting virial masses are reliable. But we must be cautious. 
For instance, there are compelling reasons to believe the BLR 
is actually a disk-wind (e.g., Murray & Chiang 1997; Proga et 
al. 2000; Proga & Kallman 2004) whose kinematics depend on 
both the mass and the accretion rate onto the BH. In this sce¬ 
nario, the virial mass calibration would depend systematically 
on BH mass and luminosity. We require objects spanning a 
wide range of AGN properties to properly calibrate the primary 
rung in our AGN BH “mass ladder.” At the same time, we may 
hope to learn about evolution of the Mbh - cr* relation by look¬ 
ing at the full distribution of local BHs (Robertson et al. 2005), 
and AGNs in particular (e.g., Yu & Lu 2004), in the Mbh-ct* 
plane. 

2. SAMPLE SELECTION AND METHODOLOGY 

Our goal is to directly compare Mbh with cr* for a large sam¬ 
ple of AGNs. As outlined in Greene & Ho (2005c), we selected 
spectroscopically identified AGNs from the Third Data Release 
of the Sloan Digital Sky Survey (SDSS DR3; Abazajian et al. 
2005) with z < 0.05 and signal-to-noise (S/N) ratios per pixel 
> 18 in the region surrounding the Ca II AA8498, 8542, 8662 
triplet (CaT). Here we consider the 32 objects from Greene & 
Ho (2005c) with resolved cr* measurements and robust Ha line 
widths. To increase the sample, we have relaxed the z require¬ 
ment and included an additional 24 galaxies to satisfy the cri¬ 
teria established by Greene & Ho for use of the “Fe region” 
(5250-5820 A): S/N > 20 (( S/N ) = 40 ±8), estimated Ed¬ 
dington ratios [Lboi/TEdd, where L E dd = 1.26 x 10 38 (M B h/M 0 ) 
ergs s -1 ] < 0.5, and AGN contamination < 80%. In addition to 
the 56 objects from the SDSS, we include 15 intermediate-mass 
BHs (IMBHs; M BH < 10 6 M 0 ) from Greene & Ho (2004) with 
cr* measurements from Barth et al. (2005), resulting in a total 
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of 71 objects (see Table 1). 

The stellar velocity dispersions are measured using a 

Table 1. Sample Properties 


Name 

(1) 

(2) 

(3*) 

FWHM Ha 

(4) 

log L Ha 
(5) 

log M bh 
(6) 

SDSS J000805.62+145023.4 

0.0454 

140 ± 27 

7610 ± 380 

41.13 ± 0.04 

7.7 ± 0.1 

SDSS J004236.86-104921.8 

0.0419 

78.4 ± 10 

1960 ± 97 

41.58 ± 0.14 

6.7 ± 0.1 

SDSS J011703.58+000027.3 

0.0456 

98.8 ± 16 

2270 ± 110 

41.45 ± 0.08 

6.8 ± 0.1 


Note. — Col. (1): Name. Col. (2): Redshift. Col. (3): cr* (km s *) measured from Ca triplet, 
unless otherwise noted. Col. (4): FWHMjj^ (km s — *). Col. (5): Ln a (ergs s — *). Col. (6): 
M B H = (2.0 +J;J) X10 6 Mq (Lhc/ 10 42 ergs g -l) 0 .55+0.02 ( FW HM Hct /10 3 km B -l)2.06±0.06 
(Greene & Ho 2005b). Note that these are formal uncertainties; actual uncertainties in Af BB are 
probably dominated by systematics in BLR geometry. Table 1 is available in its entirety via the link 
to the machine-readable version above. A portion is shown here for guidance regarding its form and 
content. 

direct pixel-fitting algorithm described in detail in Greene & 
Ho (2005c). The uncertainties are derived from simulations in 
which the AGN contamination and S/N ratio are varied for a 
grid of model galaxy spectra. 

Virial masses are based on reverberation mapping (Bland- 
ford & McKee 1982), which uses the measured lag between 
variability in the photoionizing continuum and emission lines 
to estimate BLR radii. Currently there are 35 reverberation- 
mapped AGNs in the literature (Peterson et al. 2004). These, in 
turn, are used to calibrate the radius-luminosity relation, RblrOC 
T 5100 0 ' 64 (Greene & Ho 2005b; see also Kaspi et al. 2005), from 
which it is possible to infer a radius from the AGN luminos¬ 
ity. Combining the radius with an estimate of the velocity dis¬ 
persion of the BLR yields a virial mass estimate for the BH, 
Mbh=/R v 2 /G, where / is a factor that accounts for the geom¬ 
etry of the BLR. We assume / = 0.75 for a spherical BLR (Net- 
zer 1990). Because our sample is selected to have strong stellar 
continua, virial BH mass estimation requires special care. Un¬ 
der these circumstances the most robust virial mass estimator is 
that based on the width and luminosity of the broad Ha emis¬ 
sion line, as advocated by Greene & Ho (2005b). The emission¬ 
line fitting is described in Greene & Ho (2004, 2005b). Briefly, 
the stellar continuum is modeled and removed using a prin¬ 
cipal component analysis designed for the SDSS data by Lei 
Hao (Hao et al. 2005). We then construct a multi-component 
Gaussian model of the [S II] AA6716, 6731 doublet, which is 
shifted and scaled to fit the narrow Ha+[N II] AA6548, 6583 
complex. The remaining (broad) Ha flux is fit with as many 
Gaussian components as needed to provide a statistically ac¬ 
ceptable fit, although we attach no physical significance to the 
individual components. As described in Greene & Ho (2005b), 
we measure both FWHMh 0 and the true second moment (ana) 
from the multi-component Gaussian fits. The uncertainties in 
the luminosities and line widths are derived using simulations, 
although we set a minimum uncertainty of ~ 5% on the line 
width. 


3. RESULTS 

We present the largest single collection of cr* measure¬ 
ments in active galaxies. The single-epoch Mbh values for 
the SDSS and Barth et al. (2005) samples are estimated us¬ 
ing FWHMn a (rather than a n a ; e.g., Peterson et al. 2004). 
In our full sample, we also include the 15 reverberation- 
mapped AGNs considered by either Onken et al. (2004) or 
Nelson et al. (2004), using the weighted-mean cr* from the 
two works. Weighted-mean virial products for this sample 
are computed using the H/3 lag and FWHM measurements 
presented in Peterson et al. (2004; the weights are calculated 
taking into account the asymmetric error bars), and BH masses 


are derived assuming / = 0.75. Finally, we include the well- 
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Fig. 1.— TheMeH - relation for active and inactive galaxies. Open cir¬ 
cles are measurements from this work using CaT, open squares are measure¬ 
ments using the 5250-5820 A Fe region (see Greene & Ho 2005c), and open 
triangles represent the Keck data from Barth et al. (2005) but with Mbh re¬ 
calculated using FWHMhq values from DR3. Literature data include those of 
Onken et al. (2004) and Nelson et al. (2004; crosses; BH masses from Peterson 
et al. 2004), NGC 4395 and POX 52 (Filippenko & Ho 2003; Peterson et al. 
2005; Barth et al. 2004; asterisks), and the primary sample of inactive galax¬ 
ies with dynamically determined values of Mbh (Tremaine et al. 2002; filled 
squares). Representative error bars are shown in the upper left for the SDSS 
and Keck measurements. The dashed line represents the fit for the Mbh — 
relation of inactive galaxies as given by Tremaine et al. (2002). The solid 
line is our best fit for the AGN sample with a fixed slope of (3 = 4.02 and 
a = 7.96 d= 0.03; the dotted lines show the measured intrinsic scatter of 0.4 
dex. We display histograms of the cr* (top) and Mbh (fight) distributions for 
the Keck (shaded), SDSS (open), literature (filled), and inactive (dash-dot line) 
samples. 


studied intermediate-mass BHs in NGC 4395 (reverberation 
mass Mbh= 3.6 x 10 5 Mq, Peterson et al. 2005; virial mass 
M bh = 1.4 x 10 4 M & , Filippenko & Ho 2003) and in POX 52 
(virial mass Mbh= 1.6 x 10 5 Mq, Barth et al. 2004). The full 
sample of 88 objects covers a range in mass of 10 5 - 10 8,5 Mq, 
and a corresponding range in cr* of 30 - 268 km s -1 (see Table 
1). Consistent with previous work. Figure 1 shows that active 
galaxies apparently follow the Mbh - cr* relation defined by in¬ 
active systems, even with a sample size increased by a factor of 
nearly 6. In detail, however, we do find significant differences 
in zeropoint and slope compared to the inactive sample, which 
we describe below. 

Assuming a log-linear form for the Mbh - a* relation, 
\og(M Bli /M & ) = a + /3log(ajao), with cr 0 = 200 km s -1 , we 
can formally quantify deviations of AGNs from the Mbh- a* 
relation of inactive systems. We use the Levenberg-Marquardt 
algorithm as implemented by mpfit to minimize y 2 , defined, 
following Tremaine et al. (2002), as 
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where x, correspond to a* ,, y, to Mbh./, and e, to the formal un¬ 
certainties in each measurement. In order to account for asym¬ 
metric uncertainties in each parameter, at each iteration e, is 
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selected to reflect the sign of v, -a —/3x,-. In the simplified case 
of symmetric errors, we recover the results of the Numerical 
Recipes routine fitexy (Press et al. 1992) as implemented in 
IDL. We begin by fixing the slope (3 to the best-fit value of 4.02 
from Tremaine et al. and investigate potential offsets in zero- 
point. We find, for our sample alone, a best-fit a = 7.96 ± 0.03, 
which corresponds to an offset of-0.17 ±0.07 from the value of 
a = 8.13 ± 0.06 of Tremaine et al. (consistent with Gebhardt et 
al. 2000b; Nelson et al. 2004; Onken et al. 2004; see summary 
in Table 2). The offset increases to-0.21 ±0.06 when the litera¬ 
ture values are included (using either the reverberation-mapped 
or virial mass for NGC 4395). Note that we are using formal 
uncertainties in Mbh. while the true uncertainties are probably 
dominated by our ignorance of the BLR geometry and the cor¬ 
responding uncertainty in how to extract the velocity dispersion 
of the virialized BLR from the observed line profile. Following 
Tremaine et al. and Gebhardt et al. (2000a), we estimate the 
intrinsic scatter as the value e (l that, when added to e y , results 
in a minimum Xr of 1. For both our sample alone and includ¬ 
ing the literature sample, we find an intrinsic scatter 1 of eo=0.4 
(dotted lines in Fig. 1). This is to be compared with the esti¬ 
mated scatter of 0.25-0.3 dex for the inactive sample (Tremaine 
et al. 2002). While there is probably intrinsic scatter in the un¬ 
derlying Mbh — <7* relation, we conservatively adopt 0.4 dex as 
the systematic uncertainty in the single-epoch Mbh measure¬ 
ments, so that we bracket the full range in possible uncertainty. 
If we repeat the fitting above increasing the uncertainties in the 
single-epoch BH masses by 0.4 dex in quadrature, the results 
for our sample alone are virtually unchanged, while we find 
a = 7.86 ±0.04 for the full sample, which corresponds to an 
offset of -0.27 ± 0.07. This offset is similar to that advocated 
by Onken et al., but it is statistically much more robust, being 
based on a final sample that is nearly 6 times larger. 

When we fit both the zeropoint and slope, we obtain a = 
7.84 ± 0.04 and (3 = 3.65 ± 0.13 (a = 7.89 ± 0.05, /? = 3.74 ± 
0.17) if we include (exclude) literature values, corresponding 
to slopes flatter than the value of (3 = 4.02 ± 0.32 for inac¬ 
tive sources by -0.37 ± 0.35 (-0.28 ± 0.36). Including the 
0.4 dex scatter to the single-epoch measurements yields an 
even flatter slope, with a = 7.79 ± 0.04 and (3 = 3.49 ± 0.18 
(a- = 7.68 ± 0.10, f3 = 2.95 ±0.32). It appears that much of 
the flattening is driven by objects with BH masses < 10 6 Mg, 
which rely on an extrapolation of the radius-luminosity relation 
and thus may be suspect. [As pointed out by Barth et al. (2005), 
these objects also may be biased by selection toward large Mbh 
at a given n,. | When we remove the Barth et al. objects, as well 
as POX 52 and NGC 4395, but include all other data, we find 
(3 = 4.19 ± 0.22 with the formal errors, or (3 = 3.85 ± 0.25 with 
the additional 0.4 dex uncertainty added. The fact that the slope 
is still shallower when the (more realistic) uncertainties are in¬ 
cluded gives us some confidence that the flattening is real. 

4. IMPLICATIONS AND CONCLUSIONS 

We have established, with statistical confidence for the first 
time, that the Mbh - <7* relation of local AGNs, while gener¬ 
ally similar to that of inactive galaxies, shows some significant 
differences. We find evidence for a lower zeropoint, shallower 
slope, and (probably) larger scatter. There are many competing 
effects, relating both to BLR physics and galaxy evolution. 


Table 2. Fit Parameters 



Sample 

(i) 

Uncertainties 

(2) 

a 

(3) 

0 

U) 


Our 

Formal 

7.96 ± 0.03 

4.02 


Full 

Formal 

7.92 ± 0.02 

4.02 


Our 

Systematic 

7.97 ± 0.06 

4.02 


Full 

Systematic 

7.86 ± 0.04 

4.02 


Our 

Formal 

7.89 ± 0.05 

3.74 ± 0.17 


Full 

Formal 

7.84 ± 0.04 

3.65 ± 0.13 

Our 

w/o IMBH 

Formal 

8.02 ± 0.08 

4.61 ± 0.37 

Full 

w/o IMBH 

Formal 

7.89 ± 0.04 

4.19 ± 0.22 


Our 

Systematic 

7.68 ± 0.10 

2.95 ± 0.32 


Full 

Systematic 

7.79 ± 0.04 

3.49 ± 0.18 

Our 

w/o IMBH 

Systematic 

7.66 ± 0.13 

2.88 ± 0.54 

Full 

w/o IMBH 

Systematic 

7.81 ± 0.05 

3.85 ± 0.25 


Note. — Col. (1): Data set considered; the sample “w/o 
IMBH” excludes the Barth et al. (2005) objects, POX 52, and 
NGC 4395. Col. (2): Adopted uncertainties. Col. (3): Zeropoint 
assuming log(MBH/A/©) = a + /31og(cr*/200 km s _1 ). Col. (4): 

Mbh — a * slope, fixed to 4.02 for the first 4 entries. 

that may contribute to the observed differences. If we posit 
that AGNs obey the relation of inactive galaxies exactly, then 
the scatter and zeropoint offset may be attributed to variations 
in the geometry of the BLR. This is the assumption made by 
Onken et al. (2004), who derive a statistical offset in the factor 
/ that scales all BH masses upwards relative to the case of a 
spherical BLR. If, however, as suggested by disk-wind models 
(e.g., Proga & Kallman 2004), the geometry of the BLR de¬ 
pends on physical properties of the AGN such as the BH mass 
and the Eddington ratio, then virial mass estimates will not scat¬ 
ter randomly about one fixed value of /. Rather, / will change 
systematically with the state of the system, and its average value 
for any given sample will depend on the range of parameter 
space spanned by the objects used to derive it. It is important to 
recognize that all work so far—including ours—has considered 
a rather limited range in Mbh and Lboi/EEdd, and so may lead 
to a biased value of /. Apart from the Greene & Ho (2004) 
objects studied by Barth et al. (2005), the practical challenge of 
detecting stellar absorption features to measure a,, inevitably 
biases the final sample toward relatively low Eddington ratios 
and BH masses. Excluding the Barth et al. objects, the sample 
summarized in Figure 1 has a median Lboi/EEdd = 6 x 1CT 2 and 
10 6 <j Mbh < 10 s Mg. While we expect that secondary param¬ 
eters (e.g., Lboi/EEdd) may ultimately help to account for the 
overall scatter in the AGN Mbh ~ <7, relation, we refrain from 
discussing this issue here because we believe that a proper anal¬ 
ysis would require a larger and more complete sample than is 
currently available. 

Apart from systematic uncertainties in virial masses, evolu¬ 
tion of the Mbh -o’* relation with cosmic time (e.g.. Shields et 
al. 2003; Treu et al. 2004; Walter et al. 2004; Peng et al. 2005) 
also imprints scatter and skewness into the local relation as a 
function of Mbh (Robertson et al. 2005). However, it is unclear 
at this stage how much scatter can be attributed to the fact that 
BHs in AGNs, in so far as they are radiating and thus accret¬ 
ing, are still gaining mass. Any differential growth between 
BHs and bulges will introduce additional scatter to the AGN 
Mbh — (7* relation if AGN accretion and star formation are not 
precisely synchronized (Ho 2005; Kim et al. 2005). If we as¬ 
cribe all the observed scatter to differential growth (taking the 
observed scatter of ~ 0.27 dex in the relation for inactive ob- 


The intrinsic scatter increases to 0.5 dex when an a , rather than FWHMh q is used to estimate Mbh> partly because of the difficulty of measuring the low-contrast 
line wings. We therefore consider only masses derived using FWHMh q in the following. We estimate that disk rotation contributes at most a small error to a * since 
we find no correlation between galaxy inclination angle and deviation from Mbh - 0 * (78% probability of no correlation with Kendall’s r), and typically only the 
inner ~ 20% of the galaxy light enters the fiber. See further arguments in Greene & Ho (2005a). 
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jects from Tremaine et al.), there can be no more than ~0.3 dex 
scatter (factor of 2) introduced by relative BH-bulge growth. 
This level of growth has only modest fuel requirements that are 
easy to sustain in nearby galaxies; a 10 6 M. BH requires only 
0.022 Mq yr -1 to double its mass in a Saltpeter time (4.5 x 10 7 
yr). 

Finally, taking our entire data set at face value, we do find 
evidence for a shallower slope than the inactive Mbh - cr* rela¬ 
tion. Whether this is the result of flattening at low mass or a 
different slope in AGNs is difficult to determine at the present 
time, in the absence of more AGNs with Mbh > 10 8 M. . In 
this regard, Wyithe (2005) argues that the slope for the inac¬ 
tive sample steepens considerably when the four smallest cr* 
values are removed, independently suggesting flattening at low 
mass. Flattening at low mass may be a result of the changing 
efficiency in AGN feedback in low-mass halos. For instance, 
the formalism of Vittorini et al. (2005) finds that a combina¬ 
tion of decreased optical depth (proportional to galaxy radius) 
and cooling times shorter than a dynamical time (small, dense 
halos) allows these BHs to grow to larger relative masses with¬ 
out feedback limitations. Alternatively, different growth modes 
between spheroidal and disk-dominated systems may result in 
a different final position on the Mbh-ct* plane. On the other 
hand, we cannot preclude the possibility that a change in slope 
in the BLR radius-luminosity relation at low luminosity is re¬ 


sponsible for the change in slope in the Mbh-<t* relation. How¬ 
ever, we note that a radius-luminosity slope closer to the theo¬ 
retically preferred value of 0.5 (e.g., Kaspi et al. 2000) would 
only increase the observed discrepancy in slope. Reverberation 
mapping of low-mass BHs would be required to address this 
issue. 

In summary, although we have increased the population of 
AGNs with cr* measurements by a factor of nearly 6, we do not 
yet have a large enough sample, or adequate coverage of param¬ 
eter space, to uniquely identify the cause or causes responsible 
for the observed differences between the Mbh - cr* relation of 
local active and inactive galaxies. Future effort should focus 
on (1) further enlarging the sample of AGNs with robust cr* 
measurements, (2) pushing the samples to the extremes of the 
mass distribution, particular above 10 s Mq, (3) extending the 
luminosity coverage to include objects with a wider range in 
Eddington ratios, and (4) better characterizing the BLR radius- 
luminosity relation over a broader range of AGN properties 
than is currently available. 

We thank L. Hao for providing her PCA code, D. Proga for 
useful discussion, and the SDSS collaboration for providing the 
extraordinary database and processing tools that made this work 
possible. We thank an anonymous referee for a prompt and pos¬ 
itive review. 
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